Entropic Entanglement Criteria for Continuous Variables 
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We derive several entanglement criteria for bipartite continuous variable quantum systems based 
on the Shannon entropy. These criteria are more sensitive than those involving only second-order 
moments, and are equivalent to well-known variance product tests in the case of Gaussian states. 
Furthermore, they involve only a pair of quadrature measurements, and will thus prove extremely 
useful in the experimental identification of entanglement. 
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Quantum entanglement is the property that differenti- 
ates quantum mechanical systems from classical ones. As 
such, the detection and characterization of quantum en- 
tanglement is one of the prominent goals in Quantum In- 
formation. In the discrete variable case, many detection 
schemes for entanglement have been proposed (see Q for 
review). In the continuous variable (CV) case, detection 
of entanglement is more challenging due to the compli- 
cated Hilbert space structure, and many tests only detect 
entanglement that appears in the second-order moments 
[1, 0, M S S 13 5 which is completely adequate for the 
case of Gaussian states. However, non-gaussian states 
and processes have been shown to not just enhance cer- 
tain quantum information protocols such as teleporta- 
tion [1, 0], but in fact be necessary for certain tasks, 
such as universal quantum computing [13, 11 1 and en- 
tanglement distillation 12, [13]. Towards the detection 
of CV entanglement in general, Shchukin and Vogel have 
derived an infinite hierarchy of conditions for positive 
partial transpose involving higher-order moments 
Although powerful, these conditions may not always be 
experimentally convenient [15]. 



Here we derive several entropic entanglement crite- 
ria for CVs. In contrast to previous work based on 
quantum- mechanical generalizations of entropy functions 

m E H [H, 

our criteria involve the Shannon en- 
tropy of probability distributions of a pair of comple- 
mentary quadrature measurements. We will show that 
these conditions detect entanglement in many states that 
any second-order test will not. A first set of inequali- 
ties is most sensitive, but is valid only for pure states. 
Inspired by previous work in discrete variables [lO], we 
use the entropic uncertainty relations for complementary 
CV observables [21] to derive a second set of inequali- 
ties. These have the distinct advantage that they can be 
extended to include mixed bipartite CV states. These in- 
equalities are more sensitive than the usual criteria based 
on second-order moments [H, [sl, 0, HI , and are equivalent 
to a well-known variance product criteria [3] in the case 
of bipartite Gaussian states. At the same time they 
are no more experimentally demanding than the widely 
adopted tests Si,!!. 

As in other C V inseparability criteria [1, 0, [HI , we con- 



sider the global operators 

r± = ri ± X2 



where 



S± = Si ± S2, 



cost/jXj + smt/jPj, Sj 



(la) 



(lb) 



and Xj and pj are the usual canonical variables satisfying 

1, 2 refers to each subsystem of 



i% , and j. 



a bipartite state. Note also that the operators and Sj 
satisfy [rj,Si] = i^^j. The entropy associated to a mea- 
surement of r is given by the Shannon entropy 



drR{r) In i^(r). 



(2) 



where R{r) is the probability distribution associated to 
the measurement of r, and similarly for H[S]. 
We will derive inseparability criteria of the form 



(3) 



where R± and are the probability distributions associ- 
ated to measurement of r± and s^, respectively, and c > 
is a real constant. Any separable state obeys inequality 
(|3|), while entangled states may not. For example, the 
left side of Eq. (|3]) vanishes for the common eigenstates 
of r_ and s+ or r+ and 5_, which correspond EPR-like 
states (note that [r+,5_] = [r_,5+] = 0). 

Let us first derive an inequality for the case of pure 
states. We will then extend our results to include mixed 
states. A separable pure state can be written in the 
form IV^i) |'^2), and has a corresponding wave function 
^(^1,^2) = V^i(^i)V^2(^2)- Using Eq. ^ and changing 
variables gives 



T / N 1 , / r+ + r_ 



V^2 



r+ — r_ 



(4) 



The probability distribution associated to the measure- 
ment of r± is given by 



R2 



r+ — r_ 



/ 



drRi{r)R2{Tr ± r±) = i?i * 



(5) 
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where Rj{r) = |?/;j(r)p, the symbol denotes convo- 



lution and K 



(+) 



(-) 



R2{—^)- Using the 



entropy power inequality [22|, |22 

exp(2i^[A * B]) > exp(2i^[A]) + exp(2i^[5]), (6) 

and also the fact that the Shannon entropy is invariant 
under reflections |22|, we have 



H[R±] > - In {exp{2H[Ri]) + exp(2i/ [i?2])} • (7) 
We arrive at an equivalent inequality for iJ[S'^]: 



H[S^] > -\n{exp{2H[Si])+exp{2H[S2])}. (8) 



Combining ([7j) and (|8]), we have 



H[R±] + H[S^] > 



■ (9) 



Eq. (O gives two inequalities that are satisfied by sep- 
arable pure states. Violation of either inequality (|9]) is 
then a sufficient condition for entanglement. 

We will now show that it is possible to arrive at a 
weaker pair of inequalities, and then extend them to in- 
clude mixed states. Using the entropic uncertainty rela- 
tion for continuous variables (j = 1,2) [2l| . 

H[Rj]^H[Sj] > luTre, (10) 

gives exp{2H[Rj] + 2H[Sj]) > (7re)^ which leads to 

H[R^] + H[S^] >i In 1 2(^e)2 + ^ ^i2Hm^2H[s,]) ]^ 

(11) 

Using relation (fTD|) again, one obtains 

H[R±]+H[S^] >iln{2(7re)2+ 

2(7re)2 cosh {2H[S2] - 2H[Si])} . 

(12) 

Since the hyperbolic cosine is lower-bounded by 1, we 
have 



H[R±]+H[S^] > ln27re. 



(13) 



The term on the right side is a state-independent lower 
bound for H[R±] + H[S^] of any separable pure state. 

To extend the inequalities (fT3|) to include mixed states, 
we use the fact that any bipartite separable state p can 
be decomposed into a convex sum of pure states 

p = ^Xk li^ik) ii^ikl ^ \i^2k) {i^2k\ , (14) 



where A/e > and = 1. The probability distribu- 

tions associated to a measurement r± is 



= ^ \kRk±- 



(15) 



Here Rk± is the probability to detect r± for each pure 
state in the decomposition ([El). The concavity of the 
Shannon entropy [23| gives 



H[R±]>Y,\kH[Rk±], 



(16) 



and likewise for H[S±]. As H[Rk±] = H[Rki * Rk2], 
inequality jH]) gives 



A similar condition holds for H[S±]. Summing these two 
inequalities, and using the fact that the left side of ([9]) is 
lower-bounded by ln27re, gives 



— ^ ln(27re) = In 27re, 



(18) 



which is identical to ([T3|). Thus, inequalities ([T3]) are 
satisfied by both pure and mixed separable states. We 
note also that one can take the supremum of the first 
inequality (fTSj) over all possible decompositions of the 
mixed state p to arrive at a stronger inequality. However, 
this is not suitable for experimental purposes. 

We also note that inequalities ([13]) can also be obtained 
using the positive partial transpose criterion [2] as we will 
now show. First, we note that the marginal distributions 
under partial transposition are: 



R± = R± and S± = 



(19) 



Noting that [r±, s±] = 2z, the relations in ([T9|) imply that 
any separable state must verify the uncertainty relation 
H[R±]^H[S±] > ln27re, which leads directly to inequal- 
ities (fT3|) . Furthermore, we see in ([19]) that the partial 
transposition interchanges the variables 5+ and s_, which 
is the key to obtain entanglement criteria based on un- 
certainty inequalities that can arise from the noncomu- 
tativity of r± and s±. 

An upper bound to the left side of (fTS]) can be ob- 
tained by considering that the Shannon entropy of a con- 
tinuous variable with variance is maximized when the 
probability distribution is Gaussian, for which i^causs = 
[il,!!^. Then, 



InV^ 



ln27recr±(5^ > H[R±] ^ H[S^] > ln27re. 



(20) 



where and 6'^ are the variances of R± and S± , respec- 
tively. This upper hmit is reached for Gaussian states, 
in which case we recover the Mancini-Giovannetti-Vitah- 
Tombesi (MGVT) product inequahty Q 



a±S^ > 1. 



(21) 



The left side of the double inequality ([2Q|) proves that 
the conditions (p!3|) are more sensitive than the variance 
product criteria, with equivalence in the case of Gaussian 
states. This is in accord with the extremality of entangled 
Gaussian states [2^. 

In the definition of the operators ([T]) we included the 
parameter 6j to account for local rotations of the quadra- 
ture measurements. To successfully employ the entropic 
criteria, it is necessary to find suitable rotated quadrature 
operators, parametrized by angles Oi and O2. In addition, 
one can optimize over local squeezing parameters, which 
can be included a posteriori [7[. So, in order to complete 
the analysis of the effect of real linear canonical transfor- 
mations over the single mode quadrature operators (that 
form the real symplectic group 6*^(2, R)) let us consider 
now the effect of local squeezing, which can be accounted 
for by redefining the rotated operators as 



airi ± a2r2 



Si Si ± S2. 

ai a2 



(22a) 



(22b) 



Substituting these operators, inequality (|9]) becomes 



H[R'^]+H[S':,]> - In { ^ 



{2H[Ri]+2H[Sj])+2\n - 



(23) 

and applying the entropic uncertainty relation (|10p re- 
sults in 



H[R'^]+H[S' ] > ln(27re), 



(24) 



where R'^ and S'^ are the probability distributions for 
measurements of the operators ([22]) . It is straightforward 
to show that both entropy inequalities ([23|) and ([24j) re- 
duce to (|9|) and (p!3|) when ai = a2, which demonstrates 
that these inequalities are invariant to equal amounts of 
local squeezing. 

Examples. Inequalities ([9|) and (fT3|) are more sensitive 
than criteria involving sums or products of variances. We 
will now illustrate this point with several examples of 
non-Gaussian states. Let us first consider a non-gaussian 
wave function of the form 



(25) 



This state is non-separable for all values of <j±. For this 
state both the Simon PPT criteria which is a nec- 
essary and sufficient condition to detect entanglement in 



a /a. 
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FIG. 1: (Color online) Pictoral representation of three insepa- 
rability criteria for the pure state ()25p as a function of a-/a^. 
The criteria called "Entropy inequality" and "Strong Entropy 
inequality" correspond respectively to Eqs.([T3} and (|9|. The 
dark and light grey regions correspond to the intervals of 
(7-/(7+ where each criteria detect entanglement and blank 
regions where they do not. Both criteria presented here are 
stronger than the Simon PPT condition [2|]. 



Gaussian states, and the MGVT criteria with Oi = O2 = 
0, detect entanglement provided a-/a-^ > ~ 1.732 
or cr_/cr+ < 1/V3 ^ 0.577. With Oi = O2 = 0, 
the inseparability criteria (fT3|) gives H[R±] + H[S^] = 
In (47re^cr+/cr±), where 7 = 0.577 ... is Euler's constant. 
Entanglement is detected provided a-/a-^ < e^^~^^ 12 ^ 
0.763 or cr_/cr+ > 2/e(i-^) ^ 1.310. Thus, the en- 
tropy criterion ([13]) is more sensitive than the Simon 
and MGVT conditions. Numerical results show that the 
ranges in which the pair of inequalities ([9]) detect entan- 
glement overlap, indicating that they always detect en- 
tanglement in the state ([25]) . A pictoral representation 
of these results is shown in FIG. [TJ 

We tested these criteria numerically for a number 
of combinations of low-order Fock states. A particu- 
larly interesting example is the state |0) = |0,0) /a/2 + 
|2,0) /2 + |0,2) /2, which is undetectable with any second- 
order criterion [25]. Entanglement is detectable using 
the stronger entropic inequality ([9]) with d\ = — 7r/4 and 
6>2 = 7r/4. We also tested "NOON" states of the form 
(|7V,0) + |0,iV))/\/2, where \N) is an AT-photon Fock 
state. None of these states is detected by the Simon PPT 
criteria [2] nor inequalities (p!3|) . The strong criteria Q 
detects entanglement up to A" = 5 with 6>i = ^2 = 
for all values of A" except A" = 2, for which 6^1 = and 
^2 = 7r/2 were used. 

To further test these criteria we generated uni- 
form random pure states [2^ of the form = 
Zln m=o ^nm and tested inequalities (|9|), (p!3|) and 

also the MGVT criteria, while varying angles Oi and 62 
in intervals of 7r/4. Results are summarized in table HI 
The strong inequality (|9|) detected more than 62% of the 
states up to I) = 7. 

Criteria (p!3|) also applies to mixed states. Let us con- 
sider a dephased cat state characterized by the parameter 



4 



# states 


D 




'^weak 


nMGVT 


6000 


2 


74.4% 


17.3 % 


9.9 % 


1600 


3 


86.3% 


0.5 % 


0.2 % 


800 


4 


84.9 % 


0% 


0% 


720 


5 


81.0 % 


0% 


0% 


120 


7 


62.5% 


0% 


0% 



TABLE I: Results for random states (see text), nstrong, 
Tiweak and ifiMGVT are the percentage of states detected by 
inequalities (|9]), (|13|) and ()2ip . respectively. In all cases, the 
angles and O2 were scanned in intervals of 7r/4. 



Simon [ij. The higher-order criteria i n [13 requires even 
more involved measurement schemes [15|. In summary, 
the evaluation of our entropic criteria requires the same 
experimental resources as those required to evaluate the 
commonly employed second-order inequalities in [1, 0, HI , 
while providing a more sensitive entanglement test. We 
thus expect that the inseparability tests presented here 
will be of great use in experimental settings. 

Financial support was provided by Brazilian agencies 
CNPq and FAPERJ. 




FIG. 2: (Color online) Violation of entanglement criteria ()13p 
for the dephased cat state ([26]). The vertical axis is the differ- 
ence of the left-hand side (LHS) and right-hand side (RHS) 
of ([13]). See text for details. 



< p < 1 , given by 

p =N{a) {|a, a) (a, a\ + |— a, —a) (—a, —a\ 

-(1 -p){\a,a) {-a, -a\ + -a) {a,a\)} , (26) 

where N{a) is a normalization constant. This state is 
separable only when p = 1, and is undetected by any sec- 
ond order criteria for any value of p. The entanglement 
criteria ([T3|) is shown in FIG. [2] as a function of a and p 
for real a. Using Oi = O2 = 0, the sum H[R-] -\-H[S-^] is 
less than ln27re, and thus entanglement is detected, for 
a large range of a and p. 

Let us briefly discuss the application of these entropic 
criteria in an experimental setting. For fixed values 

01 and O2 of the local rotations the Shannon entropies 
H[R±] and H[S^] can be calculated using the marginal 
probability distributions R± and S^. These can be de- 
termined directly via measurement of r± and s^, or cal- 
culated from the joint probability distributions i^(ri, r2) 
and S{si^S2). We stress that the sole determination of 
these probability distributions does not allow one to cal- 
culate arbitrary moments involving products of the rj 
and Sj operators. Hence, these measurements alone are 
not enough to determine all the second-order moments 
required to evaluate even the second-order criterion of 
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